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LONG-TERM  GOALS 

Our  research  vision  is  to  develop  and  transform  ocean  modeling  and  data  assimilation  to  quantify 
regional  ocean  dynamics  on  multiple  scales.  Our  group  creates  and  utilizes  new  models  and  methods 
for  multiscale  modeling,  uncertainty  quantification,  data  assimilation  and  the  guidance  of  autonomous 
vehicles.  We  then  apply  these  advances  to  better  understand  physical,  acoustical  and  biological 
interactions.  We  seek  both  fundamental  and  applied  contributions  to  build  knowledge  and  benefit  naval 
operations. 

A  main  focus  of  this  research  is  the  role  of  stochastic  forcing  on  ocean  uncertainty  and  variability 
predictions.  The  work  includes  collaborations  with  NRL-Stennis  to  prepare  the  transfer  of  a  subset  of 
the  capabilities  and  software  developed  by  our  Multidisciplinary  Simulation,  Estimation,  and 
Assimilation  Systems  (MSEAS)  group.  The  research  thrusts  of  interest  to  both  NRL  and  MIT,  as  well 
as  the  specific  goals  of  the  work,  are  below. 

Research  Thrusts 

Stochastic  forcing  and  uncertainty/variability  predictions 

Sensitivity  analysis  for  forecast  quality  control,  data-model  comparisons  and  data  error  models 
Multiscale  covariance  modeling  and  mapping 

Ensemble  initialization  and  generation,  towards  non-Gaussian  ensemble  initialization 

OBJECTIVES 

Our  specific  objectives  are  to: 

i)  Develop,  demonstrate  and  transfer  techniques  for  stochastic  error  modeling  and  stochastic 
boundary  forcing  for  improved  ensemble  uncertainty  predictions  with  NCOM  and  COAMPS 
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ii)  Develop  and  transfer  software  for  ocean  data  management,  quality  control  and  automated  robust 
distribution,  including  data  error-models  and  data-model  comparison  codes 

iii)  Demonstrate  and  transfer  techniques  for  multiscale  covariance  modeling  and  level-set-based 
objective  analysis  codes  for  mapping  data  in  complex  coastal/archipelago  domains 

iv)  Develop  and  demonstrate  ensemble  initialization  and  generation  schemes,  towards  non-Gaussian 
ensemble  initialization 

v)  Apply  the  above  advances  in  collaborative  sea  exercises  of  opportunity 

vi)  Strengthen  existing  and  initiate  new  collaborations  with  NRL,  using  and  leveraging  the  MIT  Naval 
Officer  education  program 

APPROACH 

The  proposed  applied  research  in  stochastic  modeling  and  ocean  uncertainty  prediction  is  linked  to  two 
growing  fundamental  fields:  prediction  and  reduction  of  uncertainties;  and,  estimation  of  properties  by 
combining  models  with  data.  From  a  fundamental  viewpoint,  uncertainty  is  characterized  by  a 
probability  density  function  (pdf).  One  of  the  aims  of  the  applied  research  and  collaborations  with 
NRL  will  be  to  improve  the  prediction  of  such  pdfs. 

WORK  COMPLETED 

Stochastic  Forcing  and  Uncertainty/Variability  Predictions 

Discontinuous  Galerkin  Finite  Element  Methods:  We  have  extended  our  3D-in-space  Discontinuous 
Galerkin  Finite  Element  (DGFE)  Python  Framework.  Within  the  Python  framework,  we  had 
previously  implemented  simple  meshing  routines,  2D  and  3D  mesh  visualization  using  Mayavi,  high- 
order  advection  schemes,  and  high-order  Hybrid  Discontinuous  Galerkin  (HDG)  diffusion  schemes. 
Now  we  have  implemented  our  novel  HDG  discretization  for  2D  or  3D  ocean  flows.  Our 
implementation  works  for  hydrostatic  or  non-hydrostatic  pressure  and  rigid-lid  or  free-surface  flows. 
For  the  free-surface  version,  we  have  mass-conserving  schemes  for  both  a  fixed  mesh  approximation 
and  a  moving  mesh.  The  theoretical  formulation  of  the  novel  HDG  scheme  has  been  improved  by 
mathematically  finding  the  correct  value  of  the  HDG  stabilization  parameter  for  the  pressure  when 
using  Projection  methods.  Our  HDG  formulation  has  also  been  extended  for  vertical  integration  and 
implicit  surface-pressure  calculation.  Additionally,  we  have  developed  and  implemented  a  new 
selective  slope  limiter  to  stabilize  the  high-order  computation,  preventing  Gibbs  oscillations.  We  have 
validated  the  advection  and  diffusion  operators,  verifying  the  correct  rates  of  convergence,  and  we  are 
currently  validating  the  ocean  model. 

Work  is  underway  on  a  parallel  DGFE  C/C++  framework,  which  is  designed  to  optimize  (and  replace) 
several  time-consuming  routines  used  in  the  high-order  HDG  scheme.  The  implicit-explicit  (IMEX) 
time  stepping  scheme  has  been  re-implemented  in  C++  and  in  parallel  using  the  PETSc  library 
package,  which  in  turn  uses  MPI  for  parallelism.  Also,  the  Numpy  C-API,  along  with  improvements 
to  the  algorithms,  has  been  used  to  speed  up  the  mesh-building  routines.  Additionally,  the  mesh  input 
procedure  to  the  model  has  been  modified  and  streamlined,  and  the  model  now  accepts  input  generated 
by  the  meshing  software  GMSH  and  distmesh. 
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Uncertain  Boundary  Conditions  and  DO  Equations:  Building  on  our  DO  equations,  we  derived  and 
implemented  a  new  efficient  scheme  for  stochastic  boundary  conditions  for  Navier  Stokes  and 
Boussinesq  equations.  The  scheme  has  now  been  validated  using  the  2D  flow  over  a  square  cylinder  in 
a  confined  channel.  We  have  conducted  a  number  of  simulations  spanning  various  Reynolds  number 
ranges,  and  examined  the  behavior  of  the  recirculation  zone  length,  and  vortex  shedding  period. 

Stochastic  Tidal  Forcing:  A  scheme  has  been  introduced  to  include  stochastic  variations  in  tidal 
forcing  for  ensemble  uncertainty  predictions  while  maintaining  the  dynamic  balances  and  conservation 
properties  of  the  original  tidal  forcing.  For  each  ensemble  member,  and  each  tidal  component,  a 
random  scale  factor  is  applied  to  the  amplitudes  of  the  tidal  elevation  and  velocity  (representing 
uncertainty  in  the  tidal  strength)  while  a  random  phase  shift  is  added  to  the  time  argument  of  the  linear 
tidal  forcing  (representing  uncertainty  in  the  timing  of  the  tides).  Using  the  same  scale  factor  and 
phase  shift  for  the  elevation  and  velocities  ensures  that  dynamic  balances  and  conservation  properties 
are  respected. 

Uncertainty  Quantification  of  Non-linear  Dynamical  Systems  for  Ocean  Modeling:  A  methodology 
and  software  have  been  developed  for  solving  a  generic  system  of  stochastic  differential  equations. 

The  goal  is  to  study  multiple  UQ  methods  and  compare  their  capabilities  and  performance  in  handling 
uncertainties  pertaining  to  ocean  modeling.  The  computations  have  been  made  efficient  using  sparsity 
of  matrices.  The  software  has  been  used  to  study  small  and  medium  dimensional  linear  as  well  as  non¬ 
linear  systems  with  stochastic  forcing  and  uncertain  initial  and  boundary  conditions.  The  main  novelty 
of  the  work  is  the  focus  on  additive  and  multiplicative  noise,  which  is  highly  relevant  to  modeling  of 
uncertainty  in  real  ocean  systems.  The  time  dependent  polynomial  chaos  method  has  been  further 
modified  to  handle  stochastic  forcing  in  full  spatial  domain  and  its  performance  in  comparison  to  the 
DO  method  has  been  analyzed. 

Ocean  Dynamics,  Modeling  and  Assimilation :  A  barotropic  velocity  feature  model  was  developed  for 
the  northern  recirculation  gyre  of  the  Gulf  Stream  near  the  continental  slope  off  the  eastern  US. 

Around  Taiwan,  a  river  discharge  model  was  developed  to  account  for  the  sudden  large  influx  of 
freshwater  from  the  island  during  and  after  the  Typhoon  Morakot,  which  occurred  in  August  2009. 

The  model  formulation  and  some  results  are  described  in  Mirabito  et  al.  (2012). 

Sensitivity  analysis  for  forecast  quality  control,  data-model  comparisons  and  data  error  models 

Uncertainty  Evaluation:  An  evaluation  was  made  of  the  real-time  uncertainty  forecasts  around  Taiwan 
made  during  the  QPE  IOP09  experiment.  Model-data  errors  were  constructed  from  the  misfits  between 
model  temperature/salinity  fields  and  SeaSoar  T/S  profiles.  Error  forecasts  for  T/S  were  made  from 
the  ensemble  standard  deviations  of  the  ESSE  forecasts.  Direct  comparisons  of  the  error  forecasts  to 
the  model-data  errors  were  made  by  visually  comparing  along-track  sections  of  both  errors,  and  by 
comparing  along-track  mean  RMS  and  bias  profiles  for  each  survey.  Probability  density  function  were 
generated  for  each  error  and  compared.  Two  skill  metrics  were  also  used  to  compare  the  pdfs: 
Kullback-Leibler  (KL)  divergence  and  logarithmic  (ignorance)  scores.  A  paper  detailing  these  results 
is  in  preparation. 

Data-model  comparisons:  Data-model  comparisons  have  been  carried  out  for  a  number  of  projects  in  a 
variety  of  regions  (Philippines,  Taiwan,  Mid-Atlantic  Bight,  etc.)  to  examine  model  sensitivities, 
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determine  appropriate  initial  and  boundary  conditions,  resolve  forcing  uncertainties  (especially 
atmospheric  forcing),  understand  implications  of  data  coverage  and  availability,  etc.  Parameters  for 
simulations  were  tuned  to  minimize  data-model  differences.  Appropriate  combinations  of  atmospheric 
flux  data,  based  on  temporal  and  spatial  resolution  and  availability,  from  different  models  were 
determined  for  the  QPE  IOP09  and  AWACS/SW06  experiments. 

Data/software  management :  A  number  of  MSEAS  software  packages  have  been  prepared  for 
distribution.  The  version  control  of  the  software  is  being  carried  out  via  SVN,  which  was  chosen  after 
an  examination  of  a  number  of  version-control  packages.  Documentation  of  the  software  is  initially 
generated  through  Doxygen  and  is  then  manually  augmented.  Software  for  ocean  data  management, 
including  multi-platform,  multi-sensor  data  quality  control  and  statistical  analysis  of  real-time 
observations,  is  being  finalized.  A  paper  detailing  this  work  (Leslie  et  al.,  2012)  is  nearing  completion. 

Multiscale  covariance  modeling  and  level-set-based  objective  analysis  codes 

The  MSEAS  Fast-Marching  Method  Objective  Analysis  (FMM-OA)  has  been  upgraded  by  the 
inclusion  of  an  option  to  use  second  order  discretizations  in  the  FMM  algorithm.  This  upgrade  also 
included  upgrades  to  the  initialization  scheme  for  the  FMM  algorithm  and  restructuring  of  the  code  for 
efficiency.  A  companion  code  was  released  to  compute  3D  streamfunction  fields  from  the  mapped 
temperature/salinity  fields.  This  code  employs  an  algorithm  to  minimize  inter-island  transports  in 
complex  archipelagos.  Additional  new  features  incorporated  into  the  FMM-OA  code,  include  options 
for  separate  data  files  for  synoptic  and  climatology  data,  explicitly  constructing  and  saving  the  large 
scale  analysis  fields  and  the  synoptic  scale  error  fields,  and  making  it  easier  for  the  user  to  provide 
their  own  structured  grid  to  the  FMM-OA.  A  manuscript  is  being  completed. 

RESULTS 

Stochastic  Forcing  and  Uncertainty /Variability  Predictions 

Discontinuous  Galerkin  Finite  Element  Methods :  As  an  example  of  results,  we  show  the  solution 
(Fig.  1)  of  the  lock-exchange  problem  (Hartel  et  al.,  2000)  which  tests  the  performance  of  solvers  with 
a  non-hydrostatic  pressure  under  density-driven  flows.  Our  solution  is  2nd  order  accurate  in  time,  3rd 
order  accurate  in  space,  uses  66  elements  in  the  vertical  and  264  elements  in  the  horizontal. 
Additionally,  our  new  selective  slope  limiter  is  active  in  this  simulation,  and  limits  oscillations  to  less 
than  0.003  %.  This  result  matches  that  of  Hartel  et  al.  (2000),  a  non-trivial  result  since  our  tests  show 
that  the  shapes  of  the  Kelvin-Helmholtz  instabilities  are  sensitive  to  spatial  resolution.  We  have  also 
finished  validating  the  advection  and  diffusion  operators,  and  their  spatial  rates  of  convergence  are 
plotted  in  Fig  2,  and  Fig  3,  respectively.  From  the  figures,  we  can  see  that  we  have  nearly  optimal 
convergence  in  all  cases.  We  are  currently  validating  the  ocean  solver,  starting  by  validating  the 
Navier-Stokes  equations. 

Uncertain  Boundary  Conditions  and  DO  Equations'.  We  tested  our  scheme  by  varying  the  inlet 
velocity  of  the  “2D  flow  over  a  square  cylinder  in  a  confined  channel”  benchmark.  The  inlet  velocity  is 
varied  by  specifying  a  uniform  probability  density  function  on  its  magnitude.  The  result  is  that  a  range 
of  Reynolds,  which  we  vary.  The  range  is  well  resolved  (equivalent  to  105  Monte-Carlo  simulations) 
and  spans  different  dynamical  regimes,  e.g.  from  a  steady  flow  with  a  recirculation  zone  to  an  unsteady 
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flow  with  vortex  shedding.  Our  DO  solver  is  capable  of  correctly  capturing  the  different  flow  regimes 
and  their  transitions,  as  can  be  seen  from  the  reconstructed  realizations  for  different  Reynolds  numbers 
(Fig.  4).  We  also  examined  the  growth  of  the  recirculation  zone  length  for  different  Reynolds  number 
ranges.  For  the  range  Re  =  [14,  26]  (Fig.  5)  and  Re  =  [28,  52]  (Fig.  6)  we  see  a  linearly  increasing 
recirculation  zone  length  of  slopes  0.0537±0.0011  and  0.0553±0.0087,  respectively.  For  the  Re  = 
[28,52]  range,  we  also  see  a  change  in  the  behavior  of  the  recirculation  zone  slope,  which  is  due  to  the 
change  in  dynamical  regime  around  Re=42.  After  Re=42,  the  flow  transitions  from  a  steady  behavior 
to  an  unsteady  vortex-shedding  behavior.  For  the  range  Re  [48,  72]  (Fig.  7)  and  Re  [70,  100]  (Fig.  8), 
the  vortex  shedding  period  decreases  with  increasing  Re.  This  decrease  is  nearly  linear,  but  clearly  has 
some  curvature. 

Stochastic  Tidal  Forcing :  The  stochastic  tidal  forcing  model  was  tested  in  data  driven  simulations 
near  Taiwan,  using  data  and  forcing  for  Aug-Sep  2009.  The  amplitude  scale  factor  was  in  the  range 
[0.8  1.2]  and  the  phase  shift  was  in  the  range  +/-  pi/8.  In  Fig.  9,  we  show  the  uncertainty  forecast  for 
the  vertically  average  velocity  in  the  upper  200m.  If  no  stochastic  variation  is  added  to  the  tidal  forcing 
(Fig.  9a)  the  velocities  in  the  Strait  of  Taiwan  and  on  the  shelf  are  strongly  constrained  and  the  forecast 
uncertainty  is  low.  With  stochastic  tidal  forcing,  the  uncertainties  in  the  shallow  regions  are  brought  in 
line  with  those  in  the  deep  regions  and  regions  of  special  sensitivity  to  tides  (northwest  tip  of  Taiwan, 
Taiwan  Bank  and  Peng-Hu  archipelago  23N,  1 19E)  are  exposed. 

Uncertainty  Quantification  of  Non-linear  Dynamical  Systems  for  Ocean  Modeling :  Uncertainty  in 
input  parameters  (including  uncertain  initial  and  boundary  conditions)  has  been  modeled  using  DO 
equations,  generalized  Polynomial  Chaos  and  Monte  Carlo  methods.  A  modification  of  the  polynomial 
chaos  method  to  handle  long  time  integration  of  non-linear  systems  with  stochastic  inputs  has  also 
been  studied.  We  found  that  the  modified  time  dependent  polynomial  chaos  method  is  able  to 
accurately  model  the  evolution  of  uncertainty  in  input  parameters  but  cannot  handle  stochastic  forcing 
in  large  systems  in  its  present  form.  We  also  found  that  since  the  time  dependent  polynomial  chaos 
method  in  its  present  form  is  defined  only  in  the  full  spatial  domain,  it  is  computationally  more 
expensive  than  DO  method,  which  can  be  implemented  on  a  reduced  time-dependent  subspace.  The 
DO  method  has  been  shown  to  be  adequate  for  modeling  both  additive  and  multiplicative  forcing  and 
will  be  used  to  study  uncertainty  quantification  in  coupled  ocean-acoustic  systems. 

Sensitivity  analysis  for  forecast  quality  control,  data-model  comparisons  and  data  error  models 

Uncertainty  Evaluation :  A  climatological  error  pdf  was  constructed  by  taking  all  the  synoptic  IOP-09 
data  prior  to  the  SeaSoar  survey  in  the  region  of  the  surveys,  computing  the  variance  as  a  function  of 
depth  and  assuming  that  the  uncertainty  is  25%  of  the  variance  and  fitting  a  normal  distribution  with 
this  uncertainty  variance.  The  uncertainty  skill  metrics  (KL  divergence  and  logarithm  score)  are  then 
computed  comparing  the  error  forecast  pdfs  to  the  model-data  error  pdfs  and  comparing  the 
climatology  pdfs  to  the  model-data  error  pdfs.  The  error  forecast  has  skill  if  its  skill  score  is  smaller 
than  the  climatology  value.  In  Fig.  10,  the  relative  KL  divergence  (KL  for  climatology  -  KL  for 
forecast)  /  (KL  for  climatology)  is  shown  for  nine  SeaSoar  survey  at  10m  depth  increments  up  to  150m 
depth.  The  temperature  error  forecast  shows  overall  skill;  during  the  last  five  surveys  (when  the  model 
was  better  tuned  to  the  2009  conditions)  the  error  forecast  almost  always  shows  50%  or  better 
reduction  in  the  KL  divergence  over  the  climatology.  The  salinity  shows  overall  skill,  but  not  as  much 
as  temperature.  Some  of  this  discrepancy  is  due  to  the  freshwater  runoff  from  typhoon  Morakot  which 
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was  not  included  in  the  real-time  simulations.  Fig.  1 1  shows  the  relative  logarithm  scores.  Since  the 
logarithm  score  computes  separate  pdfs  for  each  profile  in  the  surveys  it  is  more  sensitive  to  temporal 
variability  in  the  pdfs  (including  scales  not  represented  in  the  model).  As  a  result  the  error  forecast 
shows  a  reduced  skill  level  than  that  of  the  KL  divergence.  However,  the  skill  of  the  pdf  forecast  is  still 
superior  to  that  of  the  climatology  pdf  forecast.  Finally,  the  forecast  errors  pdfs  do  capture  temporal 
variability  from  sub-mesoscales  down  to  some  internal  wave  scale  (Fig.  12). 

Data-model  comparisons:  Data-model  comparisons  carried  out  for  the  QPE  IOP09  experiment 
indicated  that  the  climatological  temperature  and  salinity  data  utilized  for  initial  and  boundary 
conditions  were  significantly  different  in  structure  and  value  from  synoptic  conditions.  Typhoon 
Morakot  generated  a  significant  pulse  of  freshwater  not  represented  by  climatology  and  also  modified 
typical  local  temperature  structures.  The  climatological  data  was  then  modified  as  necessary  to 
provide  more  accurate  initial  and  boundary  conditions.  Data-model  comparison  also  indicated  where 
lack  of  real-time  synoptic  data  proved  critical  in  model  evolution  and  provided  an  impetus  to  acquire 
additional  synoptic  data  not  available  in  real-time.  Available  atmospheric  flux  data  sets  for  QPE 
IOP09  were  carefully  compared  and  evaluated:  flux  by  flux  and  field  by  field.  We  subsequently 
selected  what  we  determined  was  the  most  appropriate  combination  of  re-analysis  products  and 
gridded  those  on  4  different  grids  at  different  resolutions.  We  then  selected  the  highest  resolution 
among  all  available  fields  and  compiled/combined/merged  them  into  our  best  estimates  of  the 
atmospheric  forcing.  A  combination  of  18km  COAMPS  and  nested  15km/5km  COAMPS  data  sets 
were  chosen  as  the  basis  for  the  final  product  set  which  was  distributed  to  all  QPE  Pis. 

For  the  AWACS/SW06  experiment,  data-model  comparisons  led  to:  an  increase  in  vertical 
discretization  using  100  levels,  with  their  distribution  optimized  to  the  thermocline  structure  ;  the  use 
of  improved  E-P  and  direct  fluxes  from  WRF  and  NOGAPS  in  the  atmospheric  forcing;  the  correction 
of  amplitudes  in  the  diurnal  tidal  components;  and,  upgraded  initial  conditions  incorporating  synoptic 
data  and  pseudo  profiles  to  bolster  the  front,  WO  A  climatology  corrected  to  match  2006  slope 
conditions,  a  revised  shelfbreak  T/S  front  feature  model,  a  Gulf  Stream  T/S  feature  model  (based  on 
synoptic  data)  and  transport  feature  models  for  the  Gulf  Stream,  slope  recirculation  gyre  and  the 
shelfbreak  front.  Improvements  to  the  reanalysis  methodology  directly  led  to  improved  fidelity  of  the 
simulations. 

Multiscale  covariance  modeling  and  level-set-based  objective  analysis  codes 

The  MSEAS  FMM-OA  employs  an  upwind  difference  scheme  to  construct  the  shortest  ocean  distance 
between  a  data  point  and  every  grid  point.  These  distances  are  then  used  in  the  correlation  functions  to 
map  the  data  onto  the  grid.  By  comparing  the  fields  constructed  in  the  open  ocean  using  the  FMM-OA 
with  fields  constructed  using  a  Euclidean  distance  OA,  we  find  that  using  our  new  second  order 
numeric  code  reduces  the  errors  in  the  OA  fields  by  a  factor  of  2-3. 

Non-Gaussian  Data  Assimilation 

We  reviewed  and  completed  our  manuscripts  on  Gaussian  Mixture  Models  -  Dynamically  Orthogonal 
(GMM-DO)  data  assimilation  (Sondergaard  and  Lermusiaux,  2012, a, b).  We  first  derived  the  efficient 
non-Gaussian  data  assimilation  scheme,  focusing  on  a  time-dependent  stochastic  subspace  that  respects 
nonlinear  dynamics  and  captures  non-Gaussian  statistics  as  it  occurs.  The  properties  and  capabilities  of 
the  resulting  GMM-DO  filter  were  assessed  and  exemplified  by  applications  to  two  dynamical 
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systems:  (1)  the  Double  Well  Diffusion  and  (2)  Sudden  Expansion  flows;  both  of  which  admit  far- 
from-Gaussian  statistics. 

IMPACT/APPLICATIONS 

Better  understanding  and  modeling  of  physical  and  interdisciplinary  regional  ocean  dynamics  are 
essential  to  multiple  applications,  including  efficient  real-time  at-sea  research  experiments,  naval 
operations  and  coastal  seas  management.  Mathematical  and  computational  methods  and  systems  are 
necessary  to  predict  and  study  ocean  dynamics.  Scientific  progress  occurs  from  the  comparison  and 
optimal  combination  of  measurements  and  models  via  data  assimilation.  Interdisciplinary  linkages 
include  the  traditional  ocean  sciences  and  atmospheric  sciences,  but  also  new  relationships  with  other 
research  disciplines  within  the  framework  of  complex  system  sciences  and  engineering. 

TRANSITIONS 

Methods,  software  and  data  sets  were  transitioned  to  other  research  groups.  They  include:  NATO 
Center  for  Maritime  Research  and  Experimentation  (CMRE),  NCSU,  TNO,  UCSD,  UMass- 
Dartmouth,  WHOI,  MIT-OE,  NRL-Stennis. 

RELATED  PROJECTS 

Without  the  present  effort,  several  of  our  other  projects  would  not  be  feasible.  Interactions  also 
occurred  with  other  research  groups.  For  ocean  modeling,  assimilation,  adaptive  sampling  and  adaptive 
modeling,  this  involved:  MIT-OE/EAPS,  OASIS  Inc.,  NURC,  NRL-Stennis,  U.  Mass  and  TNO  (F-P. 
Lam). 
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Velocity 


Fig  1.  Density  contours  over  velocity  magnitude  for  the  Lock  Exchange  problem  with  Gr  =  1.25e6  at 
non-dimensional  time  T  =  10.  Our  solution  matches  that  of  Hartel  et  al.  (2000).  This  is  a  3rd  order 
accurate  in  space  simulation  using  66  elements  in  the  vertical  and  264  elements  in  the  horizontal. 


Fig.  2:  Convergence  of  advection  operator  for  different  orders  of  spatial  approximation.  Dashed 

lines  give  the  theoretical  convergence 


Fig.  3:  Convergence  of  diffusion  operator  for  different  orders  of  spatial  approximation.  Dashed 

lines  give  the  theoretical  convergence. 
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Fig.  4:  DO  simulation  for  the  flow  past  a  square  cylinder  in  a  confined  channel  with  uncertain 
boundary  conditions  with  Reynolds  number  range  RE=[28,  52]  (left).  Realizations  with  different 
Reynolds  numbers  reconstructed  from  the  DO  calculation.  Both  the  steady  recirculation  zone  and 
unsteady  vortex  shedding  regimes  are  captured.  In  the  steady  regime,  the  growth  of  the  recirculation 
zone  size  is  captured  and  in  the  unsteady  regime  the  change  in  vortex  frequency  is  captured. 


Recirculation  Zone  length  vs.  Reynolds  Number 
at  lime  65 
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Fig.  5:  The  non-dimensional  recirculation  zone  length  versus  Reynolds  number  for  the  simulation 

with  Reynolds  numbers  RE=[8,  32] 


10 


Recirculation  Zone  length  vs  Reynolds  Number 
at  time  65 


Fig.  6:  As  in  Fig.  5  but  with  RE=[28,  52] 


Vortex  Shedding  Period  vs  Reynolds  Number 
at  time  65 


Fig.  7:  The  vortex  shedding  period  versus  Reynolds  number  for  the  simulations  with  Reynolds 

numbers  RE=[48,  72] 


Vortex  Shedding  Period  vs  Reynolds  Number 
at  time  65 


Fig.  8:  As  in  Fig.  7  but  with  RE=[70, 130] 
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(a)  (b) 

Fig  9.  Effects  of  stochastic  tidal  forcing.  Contours  of  the  magnitude  of  the  forecast 
uncertainty  (stdv  of  the  ESSE  ensemble  spread)  for  the  vertically  averaged  velocity  in  the 
upper  200m  for  two  different  ensemble  forecasts  using  tidal  forcing,  i.e.  ESSE  ensemble  stdv 
with:  (a)  Deterministic  tidal  forcing  (b)  Stochastic  tidal  forcing.  Note  that  the  (a)  deterministic 
strong  tidal  forcing  constrains  the  velocities  on  the  shelf  and  reduces  the  forecast  uncertainty 
when  compared  to  the  ensemble  with  no  tides  (not  shown).  However,  the  ensemble  with 
stochastic  tides  maintains  and  enhances  the  uncertainties  in  shallow  regions 

( where  tides  are  the  strongest). 
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Fig  10.  Relative  KL  divergence  skill  metrics.  KL  divergence  comparing  pdfs  of  error  forecasts 
and  climatology  error  to  model-data  errors  for  (a)  temperature  and  (b)  salinity  from  nine 
SeaSoar  surveys  during  QPE IOP-09.  Relative  KL  shows  fractional  reduction  in  KL 
divergence  of  forecast  over  climatology;  positive  values  indicate  skill  in  forecast. 
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Fig.  11.  As  Fig.  10,  but  for  the  Ignorance  Score  (Logarithm  Score):  relative  logarithm 

score  skill  metrics. 
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Fig  12.  Contours  of  the  pdfs  for  salinity  error  at  40m  during  SeaSoar  survey  5  (OOOOZ-IOOOZ  Aug 
31,  2009).  The  x-axis  shows  the  profile  number,  the  y-axis  the  value  of  the  salinity  error.  The  top 
panel  shows  the  contours  of  the  error  forecast  pdfs,  the  bottom  panel  shows  the  climatology  pdf. 
The  model-data  errors  for  each  profile  are  shown  by  open  circles.  The  individual  ensemble  forecast 
deviations  by  black  dots.  The  forecast  error  pdf  show  temporal  variability  from  sub-mesoscale 

to  internal  wave  scale. 
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